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Abstract. We present the behavior of the magnetic field and activity indicators of the 

single late-type giant f3 Ceti in the period June 19, 2010 - December 14, 2010. We used 

spectropolarimetric data obtained with two telescopes — the NARVAL spectropolarimeter 

^0 ' at Telescope Bernard Lyot, Pic du Midi, France and the ESPaDOnS spectropolarimeter at 

* \ CFHT, Hawaii. The data were processed using the method of Least Square Deconvolution 

f^ which enables to derive the mean photospheric profiles of Stokes I and V parameters. We 

Oh measured the surface-averaged longitudinal magnetic field Bj, which varies in the interval 

' 0.1 - 8.2 G, the line activity indicators Call K, Ha, Call IR (854.2 nm) and radial velocity. 

By analyzing the B; variations, was identified a possible rotational period P = 118 days. 

r*j • A single, large magnetic spot which dominates the field topology is a possible explanation 

ryj ■ for the B; and activity indicator variations of /3 Ceti. 

Cu ' Key words: star: /3 Ceti - magnetic field, activity, late-type giant 

1 Introduction 
> ■ 

/3 Ceti (HD 4128, HR 188, HIP 3419) is a single late-type giant star of spectral 
qv ■ class K0 III with V = 2.04 mag and B — V = 1.02. It has started as an A-type 
T-|- \ star on the main sequence (Maggio et al., 1998). The fundamental parameters 
•^ . of (3 Ceti and their references are presented in TableU 

There is discussion in the literature about its evolutionary status: is it 
a first-ascent red giant or a clump giant star? According to Schroder et al. 
(1998) clump giants are not visible in the RASS (ROSAT All-Sky Survey). 
In this study they use data for /3 Ceti from RASS detections. On their HR 
diagram they pointed out the giant "clump region" (FigQj). It is obvious from 
this figure that the evolutionary track of /3 Ceti is far away from this region 
and the star is rather possibly situated at the base of the red giant branch. 

According to other authors f3 Ceti is a He-burning clump giant because of 
its photospheric abundances. Its low Li abundance loge(Li) = 0.01 and low 
values of the ratios C/N = 1.38 and C 12 /C 13 = 19 ± 2 (Luck & Challener, 
1995; Tomkin et al., 1975) suggest this classification. The FUV spectral prop- 
erties of /3 Ceti are very similar to other active clump stars like the Hyades 



Based on observations obtained at the Bernard Lyot Telescope (TBL, Pic du Midi, 
France) of the Midi-Pyrenees Observatory which is operated by the Institut National 
des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, 
and at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National 
Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre 
National de la Recherche Scientifique of France, and the University of Hawaii. 
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Fig. 1. HR diagram by Schroder et al. (1998). Plus symbols are RASS de- 
tections. j3 Ceti is pointed out as "+" and letter "e" with mass of 3.2 Mq. 
The giant clump region is indicated by the rectangle. 



giants 6 1 Tau and 7 Tau, and the Capella G8 III primary (Ayres et al., 1983; 
Ayres et al., 1998). Another survey of photospheric abundances of (3 Ceti was 
carried out by Sagesser &; Jordan (2005). They compared the results from 
different authors about the abundances of Fe, N, O, Mg and Si. According 
to these results, their suggestion is that this star has already undergone first 
dredge- up. 

(3 Ceti is classified as a single giant star with the highest X-ray luminosity 
logL-r = 30.2 erg/s in the solar neighborhood (d < 30 pc), calculated with a 
distance of d = 29 A pc (Maggio et al., 1998; Hunsch et al, 1996). With this 
high activity, the star is compared to Capella (a Aur), a long period RS CVn- 
type active binary, and 9 1 Tau (K0 III). A model of the atmosphere of /3 Ceti 
was performed by Eriksson et al. (1983) and they suggested the existence of 
coronal loops. Ayres et al. (2001) reported about a series of striking coronal 
flare events observed with EUVE during a period of 34 days starting on 1 
Aug 2000. 

(3 Ceti has V — R = 0.72 which places it to the left of the coronal dividing 
line in the HRD. This nearly vertical dividing line was proposed by Linsky & 
Haisch (1979) near V — R = 0.80 and it separates the stars into two groups 
- the solar-type group (V — R < 0.80) which indicates the existence of 
a chromosphere, a transition region and a corona, while the non-solar-type 
group (V — R > 0.80) shows only chromosphere lines. The existence of the 
dividing line is confirmed by the study of Simon et al. (1982) with a larger 
sample of stars. 

Using spectropolarimetric data, one can detect the presence of a magnetic 
field in plasma by the Zeeman effect (Zeeman, 1897) and measuring the 
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Table 1. The fundamental parameters of /3 Ceti and their references 



Parameter 

Mass 



Radius 



Interval of vafues 
2.8 - 3.2 M Q 



~ 15 R Q 

33.86 ± 0.16 mas 
29.38l°;g pc 

4750 - 5000 K 



2.45 - 3.05 

-0.09 - 0.13 
2.5 - 4 km/s 
60° 



References 

Maggio et al., 1998; Schroder et 
al., 1998; Gondoin, 1999; Allende 
Prieto & Lambert, 1999; Berio et 
al., 2011 

Jordan & Montesinos, 1991; Al- 
lende Prieto & Lambert., 1999; 
Sagesser & Jordan, 2005 
from Hipparcos Catalogues - Per- 
ryman et al., 1997; van Leeuwen, 
2007 

Jenkins, 1952; Ayres et al., 1998; 
Gondoin, 1999; Luck et al., 1995; 
Eriksson et al., 1983; Ottmann 
et al., 1998; Gratton & Ortolani, 
1986 

Jenkins, 1952; Jordan & Mon- 
tesinos, 1991; Eriksson et al., 
1983; Luck et al., 1995; Ottmann 
et al., 1998; Gratton & Ortolani, 
1986 

Luck et al., 1995; Gondoin, 1999; 
Alves, 2000; Ottmann et al., 1998; 
Gratton & Ortolani, 1986 
Fekel, 1997; Carney et al., 2008; 
Smith & Dominy, 1979; Gray, 
1982; Ottmann et al., 1998 
Sanz-Forcada et al., 2002, Gray, 
1989 

Jordan & Montesinos, 1991 
C. Jordan, private comunication 



Parallax and Distance 



Effective Temperature 



log g 



[Fe/H] 



Inclination angle i 
Rotational period 



199 days 
171 days 



4 Tsvetkova et al. 

Stokes parameters [I, Q, U, V]. For stars with weak magnetic fields it is 
possible to detect magnetic features only in Stokes V (circular polarization) 
because the signatures in Stokes Q and U (linear polarization) are much 
weaker than these in Stokes V. In our study for /? Ceti, we use Stokes I 
and V parameters which give information about the strength of the surface- 
averaged longitudinal magnetic field B^. 

2 Observations and data reduction 

Observations were carried out at two telescopes — the 2-m Telescope Bernard 
Lyot (TBL) at Pic du Midi Observatory, France and the 3.6-m Canada- 
France-Hawaii Telescope (CFHT) at Hawaii. 

TBL and CFHT use the twin new-generation fiber-fed echelle spectropo- 
larimeters NARVAL (Auriere, 2003) and ESPaDOnS (Donati et al, 2006) 
respectively, which in polarimetric mode have a spectral resolution of about 
65 000 and a spectrum coverage from near-ultraviolet 370 nm to near-infrared 
1050 nm in a single exposure. Stokes I (unpolarised light) and Stokes V (cir- 
cular polarization) parameters were obtained by four sub-exposures between 
which the retarders - Fresnel rhombs - were rotated in order to exchange 
the beams in the instrument and to reduce spurious polarization signatures 
(Semel et al., 1993). 

The raw data were processed using the automatic reduction software Li- 
breEsprit, developed for NARVAL and ESPaDOnS. Additional details re- 
garding the observing procedure and data reduction can be found in Donati 
et al. (1997). 

For further processing of the data we used the Least Squares Deconvolu- 
tion (LSD) technique (Donati et al., 1997). This technique enables averaging 
of several thousand absorption lines taken from one spectrum, which increases 
the signal-to-noise ratio (S/N) and thus allows us to detect weak magnetic 
signatures which would not be visible in individual lines. Using the LSD tech- 
nique, one can derive mean photospheric profiles of Stokes I and Stokes V 
from a complete echelle spectrum. 

FigE] shows the result after applying the LSD method to the spectrum of 
j3 Ceti obtained on November 20, 2010. From top to bottom it depicts: the 
circularly polarized mean Stokes V profile; a diagnostic null spectrum; the 
unpolarized mean Stokes I profile. All profiles are normalized with respect to 
the continuum intensity. The diagnostic null spectrum serves to diagnose the 
presence of spurious contributions to the Stokes V spectrum. 

The dataset includes observations in the period from June 19, 2010 to 
December 14, 2010 in which we have collected 25 spectra. First, spectra were 
automatically extracted using the LibreEsprit software and after that for the 
Zeeman analysis we used the LSD method with a mask calculated for an 
effective temperature T e ff = 5000 K, logg = 3.0 and a microturbulence of 2 
km/s, consistent with the literature data for the star (Thevenin, 1998; Hekker 
& Melendez, 2007). For each spectrum we calculated mean Stokes I and V 
profiles averaging about 12 700 lines. Then we computed the longitudinal 
magnetic field B; in G, using the first-order moment method (Donati et al., 
1997; Rees & Semel, 1979; Wade et al., 2000). We also measured the radial 
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LS deconvolved Zeeman signature (12720 lines used) 




-200 



-100 



Velocity (km/s) 



Fig. 2. LSD profiles of /3 Ceti observed on November 20, 2010 with ES- 
PaDOnS. From top to bottom: mean Stokes V profile; diagnostic null spec- 
trum; mean Stokes I profile. See details about these profiles in the text. 
Profiles are shifted vertically; the Stokes V profile and null spectrum are 
expanded by a factor of 200 for clarity. 



velocity (RV) on the LSD Stokes I profiles using a Gaussian fit. The RV 
stability of NARVAL and ESPaDOnS is about 30 m/s (Moutou et al., 2007). 
For better understanding of the stellar activity, we measured the activity 
indicators Ha at 656.3 nm and the Call IR triplet at 854.2 nm (relative 
intensities with respect to the continuum) and Call K (the relative intensity 
of the emission core relative to the intensity at A = 355 nm, I/I(395nm)). 

In order to determine the rotational period of the star, we applied to our 
dataset a Lomb-Scargle period search analysis (Lomb, N.R., 1976; Scarglc, 
J.D., 1982) for B; and some of the activity indicators. 

3 Results and Discussion 



The spectra allow us to obtain simultaneously the longitudinal magnetic field 
B;, the activity indicators and the radial velocity so that we can study the 
correlations between them. The list of observations and measured values of 
B; in gauss, Ha, Call K, Call IR and their uncertainties, as well as the radial 
velocity measurements, are listed in Tabl^2j 

The changes of the B/ values are in the interval from 0.1 G to 8.2 G and 
B; remains of positive polarity for all our observations. The time variation 
of B/ shows two types of behavior (Figj3]). A sine wave is clearly visible in 
the first part of our observations in the period June 19, 2010 - October 20, 
2010. We observe a plateau after the sinusoidal curve, which covers the period 
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November 12, 2010 - December 14, 2010. The existence of this second part of 
the dataset, which is completely different from the first one, is an evidence 
for a change in the surface magnetic field configuration of the star. 

The behavior of the activity indicators and their uncertainties are pre- 
sented in Fig2] and Tabk[2j Their variations are in the intervals: 0.20 - 0.22 
for Ha; 0.17 - 0.23 for Call K; 0.13 - 0.16 for Call IR. The activity indicators 
Call K and Call IR vary in nearly the same way but Ha shows no significant 
variations with time, taking into account the accuracy of its measurements. 
Comparing the curves for B; and the activity indicators Call K and Call IR 
in FigJH] and FigH]but only for the first part of the dataset, one can notice 
that the variability shows some similarity. This sinusoidal curve for B/ could 
be due to a single large spot at the stellar surface but a more complex topol- 
ogy is also possible. Moreover Call K and Call IR also vary in a manner 
similar to B/ which could support the idea of one big spot at the stellar sur- 
face which dominates the behavior of all magnetic activity indicators. It is 
not the same case for the second part of the dataset where Ha, Call K, Call 
IR show different behaviors regarding B^. For this second part, the intensity 
of the indicators decrease which means that we detect less emission and the 
chromosphere is less heated. A smaller part of the hemisphere is occupied by 
the magnetic field. 

Radial velocity variations with time are shown in Figj5j They vary in the 
interval 13.24 - 13.47 km/s. 
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Fig. 3. Time variations of the surface-averaged longitudinal magnetic field 
B/ and its uncertainty in G. 



We applied a Lomb-Scargle method for the period search of the data. 
Because of the two different behaviors of the dataset, it was not possible to 
perform the periodogram analysis with all the data. So we used only the first 
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Table 2. Data for B;, activity indicators and their accuracy 



Date 


HJD 


Call K 


a 


Ha 


a 


Call IR 


a 


Hi, [Gj 


a, |G| 


19 Jun 10 


2455368.123 


0.191 


0.005 


0.211 


0.002 


0.134 


0.003 


4.4 


0.4 


21 Jun 10 


2455370.134 


0.194 


0.004 


0.214 


0.001 


0.135 


0.003 


3.9 


0.4 


16 Jul 10 


2455395.140 


0.210 


0.004 


0.220 


0.001 


0.141 


0.001 


7.4 


0.3 


17 Jul 10 


2455396.132 


0.216 


0.005 


0.217 


0.001 


0.144 


0.002 


8.1 


0.4 


25 Jul 10 


2455404.138 


0.207 


0.012 


0.215 


0.001 


0.146 


0.000 


8.2 


0.9 


04 Aug 10 


2455414.006 


0.218 


0.005 


0.213 


0.001 


0.147 


0.001 


7.2 


0.4 


06 Aug 10 


2455415.633 


0.197 


0.007 


0.211 


0.002 


0.149 


0.001 


7.0 


0.5 


16 Aug 10 


2455425.685 


0.201 


0.004 


0.221 


0.002 


0.148 


0.003 


4.5 


0.5 


02 Sep 10 


2455442.699 


0.189 


0.007 


0.218 


0.001 


0.147 


0.001 


2.6 


0.4 


19 Sep 10 


2455459.473 


0.187 


0.005 


0.211 


0.002 


0.142 


0.001 


1.2 


0.4 


26 Sep 10 


2455466.500 


0.198 


0.006 


0.215 


0.001 


0.147 


0.002 


0.1 


0.4 


05 Oct 10 


2455475.518 


0.212 


0.007 


0.215 


0.002 


0.155 


0.002 


1.8 


0.7 


13 Oct 10 


2455483.475 


0.230 


0.005 


0.213 


0.001 


0.165 


0.001 


3.6 


0.4 


15 Oct 10 


2455485.832 


0.216 


0.003 


0.223 


0.002 


0.155 


0.002 


3.9 


0.4 


17 Oct 10 


2455486.863 


0.222 


0.006 


0.217 


0.001 


0.154 


0.001 


3.8 


0.3 


20 Oct 10 


2455490.471 


0.211 


0.016 


0.213 


0.001 


0.158 


0.001 


5.1 


0.8 


12 Nov 10 


2455513.360 


0.199 


0.004 


0.212 


0.002 


0.149 


0.004 


4.9 


0.7 


15 Nov 10 


2455516.831 


0.202 


0.003 


0.222 


0.002 


0.142 


0.003 


4.0 


0.4 


20 Nov 10 


2455521.798 


0.197 


0.007 


0.212 


0.002 


0.137 


0.001 


4.2 


0.5 


21 Nov 10 


2455522.718 


0.195 


0.006 


0.212 


0.004 


0.137 


0.004 


3.9 


0.4 


26 Nov 10 


2455527.368 


0.185 


0.004 


0.212 


0.001 


0.136 


0.003 


4.5 


0.4 


27 Nov 10 


2455528.855 


0.196 


0.006 


0.212 


0.002 


0.137 


0.002 


3.9 


0.4 


04 Dec 10 


2455535.349 


0.175 


0.009 


0.207 


0.001 


0.138 


0.003 


4.8 


0.6 


12 dec 10 


2455543.383 


0.175 


0.010 


0.211 


0.001 


0.139 


0.002 


5.0 


0.8 


14 Dec 10 


2455545.292 


0.167 


0.013 


0.208 


0.002 


0.147 


0.003 


5.3 


1.1 



Note: Relative intensity values are given for Ha and Call IR. For Call K emission core 
1/1(395 nm) is measured. B; values and their uncertainties are given in gauss. 




Fig. 4. Time variations of the activity indicators Ha, Call K and Call IR 
from top to bottom, respectively, and their uncertainties. 
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Fig. 5. Time variations of the radial velocity. The stability of the instruments 
is about 30 m/s (Moutou et al., 2007). 
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Fig. 6. Results from the periodogram analysis for Bj {left panel) and for 
Call K and the radial velocity {right panel) after applying the Lomb-Scargle 
method. 
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Fig. 7. Variations of B; with respect to the phase with the determined pos- 
sible period P = 118 days. Observations are in the period June 19, 2010 - 
October 20, 2010. Error bars are indicated. Black squares indicate the first 
rotation of the star, blue dots indicate the second rotation of the star. 
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Fig. 8. Variations of Call K with respect to the phase with the determined 
possible period P = 85 days. Observations are in the period June 19, 2010 - 
October 20, 2010. Error bars are indicated. Black squares indicate the first 
rotation of the star, blue dots indicate the second rotation of the star. 
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Fig. 9. Variations of the radial velocity with respect to the phase with the 
determined possible period P = 78 days. The stability of the instruments is 
about 30 m/s (Moutou et al., 2007). Observations are in the period June 19, 
2010 - October 20, 2010. Black squares indicate the first rotation of the star, 
blue dots indicate the second rotation of the star. 



part of the dataset. We found possible periods for B^, Call K and the radial 
velocity. 

We determined the following possible periods: for B; a period of P = 118 
days with false alarm probability (FAP) of 5%; for Call K a period of P = 
85 days with FAP of 51% and for the radial velocity a period of P = 78 
days with FAP of 18%. These periods and FAP values are shown in FigEl 
However, one should be cautious about these periods because our dataset 
covers less than two periods. 

As it was mentioned, B; is a surface-averaged longitudinal magnetic field 
and it originates from the star's photosphere. The activity indicators Ha, 
Call K and Call IR originate from the star's chromosphere. Our period 
determinations may suggest that the chromosphere of j3 Ceti rotates faster 
than the photosphere. Another interesting result is the RV period P = 78 
days which also reflects the photospheric rotation and it is shorter than the 
B; period. We will continue to obtain data for this star in order to be able 
to understand better these differences. 

Using these periods, we made new plots with the behavior of B/, Call K 
and RV with respect to the phase — Fig|7J Fig|8]and FigEJ The X-axis values 
for these three figures were calculated according to the determined periods 
respectively for B;, Call K and RV. These periods fit well to our observational 
data especially for Call K and RV. But the data do not cover completely 
at least two rotations so the periods we found are possible periods. More 
observational data and their analysis could confirm or reject these periods. 
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Conclusions 

We collected observational data for 6 months (June 19, 2010 — December 14, 
2010) for the single late-type giant f3 Ceti. The Least Squares Deconvolution 
method was applied for magnetic field detection and measurements. Our 
results are: 

- B; remains of positive polarity for all our observations with variations 
between 0.1 G and 8.2 G. 

- We found that during the last month of our observational period a change 
in the configuration of the magnetic field at the stellar surface might 
occurs. 

- The behavior of the activity indicators follows the behavior of B; but 
only for the first part of our dataset. The case is not the same for the 
second part, where the emission decrease as a result that the magnetic 
field covers a smaller fraction of the stellar surfece. 

- Using the Lomb-Scargle method for the period search, we found possible 
periods for B^ (P = 118 days), Call K (P = 85 days) and for the radial 
velocity (P = 78 days). 

Our next steps for studying (3 Ceti are: 

- to compute a Zeeman-Doppler Imaging map of the star 

- to collect more observational data which will help us to analyze further 
the magnetic variability of (3 Ceti. 
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